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ABSTRACT 

Recent numerical simulations of the fragmentation of primordial molecular clouds in hierarchical cos- 
mogonies have suggested that the very first stars (the so-called Population III) may have been rather 
massive. Here we point out that a numerous population of massive black holes (MBHs) - with masses 
intermediate between those of stellar and supermassive holes - may be the endproduct of such an episode 
of pregalactic star formation. If only one MBH with m, > 150 M Q formed in each of the 'minihalos' 
collapsing at z ~ 20 from 3-cr fluctuations, then the mass density of Pop III MBHs would be comparable 
to that of the supermassive variety observed in the nuclei of galaxies. Since they form in high-er rare 
density peaks, relic MBHs are predicted to cluster in the bulges of present-day galaxies as they become 
incorporated through a series of mergers into larger and larger systems. Dynamical friction would cause 
> 50(m./150 M0) 1 / 2 such objects to sink towards the center. The presence of a small cluster of MBHs in 
galaxy nuclei may have several interesting consequences associated with tidal captures of ordinary stars 
(likely followed by disruption), MBH capture by the central supermassive black hole, gravitational wave 
radiation from such coalescences. Accreting pregalactic MBHs may be detectable as ultra-luminous, 
off-nuclear X-ray sources. 

Subject headings: black hole physics - cosmology: theory - early universe - galaxies: formation 



1. INTRODUCTION 

The very first generation of stars must have formed out 
of unmagnetized, metal-free gas, characteristics leading 
to a well-defined set of cosmic initial conditions, a signif- 
icant simplification of the relevant physics, and an initial 
stellar mass function (IMF) which may have been very dif- 
ferent from the present-day Galactic case (Larson 1998). 
The primordial IMF is expected to play a crucial role in 
determining both the impact of Pop III stars on the ion- 
ization, thermal, and chemical enrichment history of the 
intergalactic medium (Couchman & Rees 1986), and the 
direct observability of the universe at the end of the 'dark 
ages', e.g., through transient sources such as supernovae 
or 7-ray bursts (Barkana & Loeb 2001; and references 
therein). In the absence of metals, it is molecular hydro- 
gen (H2) which acts as the main gas coolant in some of 
the earliest non-linear cosmic structures expected in hier- 
archical clustering theories, subgalactic fragments ('mini- 
halos') with total masses M > a fewxl0 5 M Q collapsing 
from high-cr fluctuations at redshift z > 20. Recent three- 
dimensional numerical simulations of the fragmentation of 
primordial molecular clouds in these popular cosmological 
scenarios - all variants of the cold dark matter (CDM) 
cosmogony - have shown convergence towards a thermal 
regime where the gas density is ~ 10 4 cm -3 and the tem- 
perature is a few hundred Kelvin (Bromm, Coppi, & Lar- 
son 1999, 2001; Abel, Bryan, & Norman 2000). Below 
this temperature the ro-vibrational line cooling time by 
H2 is longer than the collapse time, while above this den- 
sity the excited states of hydrogen molecules are in local 
thermodynamic equilibrium and the cooling time is nearly 
independent of density. At this point collapse becomes 



quasistatic rather than a runaway free-fall, and all simu- 
lations show the formation of one or several Jeans unstable 
clumps with masses exceeding a few hundred solar masses. 
Because of the slow subsonic contraction, the gas distri- 
bution is found to remain smooth prior to the onset of 
three body H2 formation (densities in excess of 10 8 cm -3 ), 
and no further fragmentation into sub-components is seen. 
This suggests that the very first stars may have been very 
massive (Bromm et al. 1999; Abel et al. 2000), as or even 
more massive than the Pistol Star, m > 200 M Q (Figer 
et al. 1998). 

In practice, the characteristic mass of the first stars and 
the fraction of the halo gas which is available for early 
star formation will depend on radiative line transfer ef- 
fects which are difficult to simulate in multi-dimensions, 
the poorly understood physics of stellar radiative and me- 
chanical 'feedback' mechanisms, and perhaps also on the 
chemistry of HD and LiH molecules (Galli & Palla 1998). 
If each clump in the above simulations ends up as a single 
very massive star (VMS) with m > 250 M Q , then a nu- 
merous population of primordial black holes with masses 
m, « m/2 > 150 M Q is the expected endproduct of such 
an episode of pregalactic star formation (e.g., Bond, Ar- 
nett, & Carr 1984; Larson 1999; Schneider et al. 2000; 
Fryer, Woosley, & Heger 2000). The early work by Bond 
et al. (1984) was motivated by the possibility that the 
remnants of VMSs could be an important contributor to 
the dark matter in galaxies (see also Carr 1994). The 
holes envisaged here, in contrast, collectively contribute a 
very small fraction of the nucleosynthetic baryons. Since 
they form in high-cr peaks, these relic, massive black holes 
(MBHs) are predicted to be strongly 'biased', i.e. more 
abundant in the central parts of present-day galaxies. The 
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properties and detectability of primordial MBHs are the 
main focus of this Letter. 

2. THE FATE OF VMSs 

The evolution and stability of zero- metallicity VMSs 
(Z = 0;m = 120-300 M©) have been recently investi- 
gated by Fryer et al. (2000) and Baraffe, Heger, & Woosley 
(2001). Above 100 M©, main sequence stars of Pop I com- 
position are known to be vibrationally unstable to nuclear- 
powered radial pulsations, an instability which leads to 
substantial mass loss. In Z = VMSs, however, the effi- 
ciency of these instabilities appears to be greatly reduced 
relative to their metal-rich counterparts due to their hot- 
ter cores. Since, at zero metallicity, mass loss through 
radiatively-driven stellar winds is also expected to be neg- 
ligible (Kudritzki 2000) , Pop III stars may actually die los- 
ing only a small fraction of their mass. If they retain their 
large mass until death, VMSs with 100 < to < 250 M Q 
(helium cores below about 150 M Q ) will encounter the 
electron-positron pair instability and disappear in a giant 
nuclear-powered explosion, leaving no compact remnants. 
In still heavier stars, however, oxygen and silicon burning 
is unable to drive an explosion, and complete collapse to 
a black hole will occur (Bond et al. 1984). 

3. FORMATION OF PREGALACTIC MBHs 

While considerable uncertainties still remain on the 
characteristic mass or mass spectrum of the first luminous 
objects, there is growing theoretical basis for considering 
the possibility that a population of metal-free stars with 
masses larger than 250 M© formed in the minihalos that 
populated the universe at high redshifts. After 2 Myr such 
a star would have collapsed to a MBH containing at least 
half of the initial stellar mass (Fryer et al. 2000). 

In a CDM flat cosmology with Om = 1, -Hio = 
100 h kms" 1 Mpc -1 , and rms mass fluctuation normal- 
ized at present to erg = 0.63 on spheres of 8/i _1 Mpc, 
minihalos of mass M = 3 x 10 5 hr 1 M© 4 would be col- 
lapsing at z = 22 from 3-a fluctuations. A halo of mass 
M at this epoch is characterized by a virial radius (i.e. 
the radius of the sphere encompassing a mean overden- 
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sity of 200) r v i r 

V c — 5 kms -1 Mg_g a at r v ; r . The gas collapsing along 
with the dark matter perturbation will be shock heated 



to the virial temperature T v - 



1500 K M 5 2/ 5 3 ^, where 
M5.5 is the halo mass in units of 3 x 10 5 /i _1 M© and /1 is 
the mean molecular weight. The total gas mass fraction 
within the virial radius is equal initially to the baryonic 
matter density parameter Clb, with Q b h 2 = 0.02. Also, in 
a Gaussian theory, halos more massive than the 3-tr peaks 
contain 0.3% of the mass of the universe. Assuming only 
one black hole of mass to. > 150 M© forms in each of the 
3-a minihalos, one estimates that a fraction 
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fraction of the supermassive variety found in the nuclei of 
most nearby galaxies, 
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of the nucleosynthetic baryons would end up in pregalac- 
tic MBHs. Interestingly, this is comparable to the mass 

4 This is well above the cosmological (gas+dark matter) Jeans mass for an Einstein-de Sitter universe, Mj 
[(1 + z)/20] 3 / 2 . It corresponds instead to the minimum mass threshold for baryonic condensation: above this mass the H2 cooling time is 
shorter than the Hubble time at virialization, the baryonic material in the central regions becomes self-gravitating, and stars can form. Recent 
numerical simulations by Fuller & Couchman (2000) find a minimum collapse mass of m 5 x 10 5 Mq (h = 0.65) at 2 = 20. 



where p C rit is the critical density today and we have taken 
for p sm bh the value recently inferred by Merritt & Fer- 
rarese (2001). If an extreme top-heavy IMF is linked, 
as suggested by numerical simulations, to primordial H2 
chemistry and cooling, then it seems unlikely that the for- 
mation of MBHs from zero-metallicity VMSs could be a 
much more efficient process, since: 

(1) Metal-free VMSs emit about 10 photons per sec- 
ond per solar mass of stars in the Lyman-Werncr (LW, 
11.2-13.6 eV) bands (Bromm, Kudritzki, & Loeb 2000). 
Over a main-sequence lifetime of 2 x 10 6 yr, 10 5 LW pho- 
tons will be radiated per stellar baryon, or approximately 
(10 5 x / m bh x 2) > a few LW photons per baryon in the 
universe from the 3-a peaks alone. Radiation in the LW 
bands will photo-dissociate H 2 elsewhere within the host 
halo via the Solomon process (Omukai & Nishi 1999), es- 
cape into the intergalactic medium (IGM), and form a soft 
UV cosmic background which suppresses molecular cool- 
ing throughout the universe (Haiman, Rees, & Loeb 1997), 
except perhaps in the cores of already collapsing clouds 
where H2 self-shielding may be important. This back- 
ground would inhibit the formation of a much more nu- 
merous population of VMSs, e.g. in M = 3 x 10 5 hr 1 M© 
minihalos at z = 14 when they collapse from 2-a fluctua- 
tions (containing 5% of the mass of the universe). 

(2) The same VMSs in the 3-a peaks also emit similar 
amounts of Lyman-continuum photons, again about 10 5 
per stellar baryon (this is > 20 times higher than the corre- 
sponding number for a stellar population with a standard 
Salpeter IMF). Photoionization by an internal UV source 
or a nearby external one will heat the gas inside a mini- 
halo to a temperature of 10 4 K, leading to the photoevap- 
oration of baryons out of their hosts. A few H-ionizing 
photons per nucleosynthetic baryon may not be enough 
to photoevaporate all minihalos and reionize the universe 
by these early times (Haiman, Abel, & Madau 2001), an 
event which must likely wait for the formation, collapse, 
and cooling of the more numerous 2 and 1-a peaks or 
of halos with higher mass. After the reionization epoch 
further fragmentation and star formation is only possible 
within more massive halos with virial temperatures T v ; r > 
a few xl0 4 K, where pressure support is reduced and gas 
can condense due to atomic line cooling. 

(3) Stars with masses in the range 100 < to < 250 M© 
are predicted to make pair-instability supernovae with ex- 
plosion (kinetic) energies of up to 10 53 erg. This is typi- 
cally much larger than the baryon binding energy of a sub- 
galactic fragment, E b w 10 50 erg hr 1 Q b M 5 5/ 5 3 [(1 + z)/20]. 
Minihalos will then be completely disrupted ('blown- 
away') by such energetic events, metal -enriched material 
will be returned to the IGM (e.g. Madau, Ferrara, & Rees 

: 10 38 h- 1 M (Q 6 h 2 /0.02)- 3 / 5 
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2001) and later collect in the cores of more massive halos 
formed by subsequent mergers, where a 'second' genera- 
tion of stars may now be able to form with an IMF which 
is less biased towards very high stellar masses. 

4. PREGALACTIC MBHs IN GALAXY BULGES 

We have seen that only a small fraction of the baryons 
need condense into VMSs in order to reheat the IGM and 
choke off the further formation of similar objects. The 
exact threshold depends on uncertain details, but the esti- 
mates in the previous section suggest that only the baryons 
within 3-cr peaks would be involved. The expected abun- 
dance of MBHs - assuming they behave similarly to col- 
lisionless dark matter particles - is then a function of the 
environment, being higher (positively 'biased') within a 
massive galaxy bulge that collapsed at early times, and 
lower ('antibiased') within a dwarf galaxy collapsing at the 
present-epoch (White & Springel 2000; Miralda-Escude 
2000). This effect can be quite easily quantified within 
the extended Press-Schechter formalism by computing the 
average number of progenitor minihalos at redshift z in a 
unit range of In M which will have merged at a later time 
zq into a more massive halo of mass Mq , 



dN 
d\nM 



D 



TT 



1J 

Mq crjvf -gy^ exp 



LP 
2S 



da 



M 



dM 



(3) 



where D — S c (z — Zq) and S = a 2 M — <j 2 Mq (Bower 1991; 
Lacey & Cole 1993). Here a 2 M (cr%[ ) is the variance of 
the linear power spectrum at z = smoothed with a 'top- 
hat' filter of mass M (Mq), and S c = 1.686 is the usual 
critical overdensity for spherical collapse in an Einstein- 
de Sitter universe. Denoting with v = 5 C (1 + z)/ctm and 
vq = 5 c (1+zq)/(Jm the 'heights' of the final density field in 
units of the rms, the halo number density can be rewritten 
as 
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where a = aM / a M < 1. A minihalo will collapse 
sooner if it lies in a region of larger-scale overdensity, 
av/2, leading to an enhanced abundance of MBHs 



> 



with respect to the mean. The bias is expected to be 
weak, e.g., in a 'Milky Way' halo (M = 10 12 /i _1 M Q , 
cm = 3) if this collapsed at zq = 0.8 from a l-er fluc- 
tuation {vq = 1). In this case the abundance of M = 
3 x 10 5 hr 1 Mq minihalos formed at z = 22 from 3-cr fluc- 
tuations (v — 3), that become part of our model Milky 
Way at zq, would actually correspond to a density peak 
equal to (y — o.vq)/(1 — a 2 ) 1 / 2 = 2.84 times the rms, or an 
enhancement of only a factor 1.6 relative to the mean. As- 
sume, on the other hand, that galaxy bulges form at early 
times. In a bulge of total mass Mo = 2 x 10 11 h^ 1 M Q 
(cmo = 3.9) collapsing at zq = 3.5 from a 2-er fluctuation, 
the abundance of minihalos would actually correspond to 
a density peak equal to 2.5 times the rms, or an enhance- 
ment factor of about 4. Figure 1 shows the number of 
minihalos - hence of MBHs - of mass M formed at differ- 
ent redshifts z which end up as part of a Mq system at zq, 
for Oo,M ) = 0.8, lO 12 /!" 1 M ('Milky Way' halo) and 
(z , M ) = 3.5, 2 x 10 11 h- 1 M Q (older 'bulge'). 

5 This is similar to the distribution of stars observed in the inner few 
& Townes 1994). 



5. PREGALACTIC MBHs IN GALAXY NUCLEI 

We have shown that, if current structure formation the- 
ories are correct and the primordial IMF was extremely 
top-heavy, a large number of MBHs - the remnants of pre- 
galactic VMSs born at very high redshifts in rare high-cr 
peaks - is expected to populate the bulges of present-day 
galaxies. As they get incorporated through hierarchical 
merging into more massive halos, these MBHs will be sur- 
rounded by dense baryonic (stars+gas) cores which will 
be gradually dragged into the center by dynamical friction 
against the visible and dark matter background. Dynami- 
cal friction will favor accretion of larger mass subunits and 
infalling satellites with low initial angular momenta (high 
eccentricities): the cores will merge and undergo violent 
relaxation. 




Fig. 1. — Mass function of minihalos formed at z = 15,20,25 
which, by the later time zq, will have merged into a more mas- 
sive halo of total mass Mq. Solid curves: zq = 0.8, Mo = 
10 12 h,- 1 Mq ('Milky Way' halo). Dashed curves: z = 3.5, 
M = 2 X 10 11 h' 1 M (older 'bulge'). 



Consider now the fate of a population of ~ 5,000 MBHs 
in the bulge of a 'Milky Way' galaxy. We assume that 
the bulge mass - Mbulge = 2 x 10 10 M© within a ra- 
dius rbuige = 1-5 kpc - is dominated by stars distributed 
with a singular isothermal profile of circular speed V c , 
p-k(r) — V 2 /(4nGr 2 ), 5 truncated at rbuigo- In our model 
this density profile is composed of stars, stellar remnants 
(including low-mass black holes), and primordial MBHs. 
Because of their large mass, MBHs 'relax' much faster than 
ordinary stars and mass segregation occurs owing to dy- 
namical friction. A hole at radius r — 10 rio pc with veloc- 
ity V c = 250 V250 kms -1 will sink towards the center on 
the timescale 
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(e.g. Binney & Tremaine 1987; Begelman, Blandford, & 
Rees 1980). Equation (|^) with a Coulomb logarithm of 
In A = 14 predicts that any MBH whose original orbit 
lays within raf = 20 (m./150 Mq) 1 / 2 pc will spiral into 

hundred pc of the Galactic center, p* cx r — 18 (e.g. Genzel, Hollenbach, 
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the center in less than a Hubble time. A central cluster 
of N mhh = 5, 000 (r df /r bulg0 ) > 50 (to./150 M ) 1 /2 MBHs 
will form with a rate of infall which decreases as i -1 / 2 and 
is currently of order two holes per Gyr per galaxy. Here the 
precise numbers depend on the stellar density profile and 
velocity distribution in the inner regions, and on the hole 
mass. Any increase in dynamical friction due to extra ma- 
terial (gas or dark matter) associated with the hole, either 
at the present-epoch or - perhaps more likely - during the 
assembly of the halo from smaller units, would enhance the 
central concentration and the number of MBHs near the 
nucleus. 

The presence of a small central cluster of pregalactic 
MBHs would have several interesting observational con- 
sequences, which will be discussed in details in a sub- 
sequent paper. Here we just want to point out that a 
close encounter may result in the tidal capture of an or- 
dinary main-sequence star of mass to* and radius i?* 
by a MBH, and the formation of a close binary sys- 
tem. Tidal capture occurs in galaxy cores whenever a 
star has a pericenter passage at r p = xR* very close to 
the tidal-breakup radius rt = (to./to*) 1 / 3 i?*. The tidal 
capture cross-section can be computed from the require- 
ment that the energy stored in the tides, « Gm 2 /(R ir x 6 ) 
is larger than the relative energy of the pair at infin- 
ity (equal to TO*i> 2 el /2 for to* <§C to,), which is only a 
small fraction of the kinetic energy available at closest 
approach. This condition yields x < (m.^/m^d) 1 / 3 , 
where v± — (2Gm i ,/ R^) 1 / 2 is the escape velocity at stel- 
lar surfacer. The cross-section leading to an encounter 
closer than x is £ = ttR 2 (t^/frci) 7 ^ 3 (mt/m*) 4 ^ 3 (Fabian, 
Pringle, & Rees 1975), in good agreement with that de- 
rived by a detailed computation of the amount of energy 
deposited into oscillatory modes during a close encounter 
with an = 1.5 polytrope (Kim & Lee 1999). For an 
isothermal stellar distribution with one-dimensional ve- 
locity dispersion a, the tidal capture rate at the galaxy 
center is then 

P*(0)iV m bh f°° r^/ \ / V rc\^j3 (r\ 
r= Ta 2W2 W re lS Kel)exp(--^)d 3 U rel . 6 

With g = V c /V2 = 180 kms" 1 , to* = 0.5 M , i?* = 
0.5 R , and /9*(0) = 10 6 M©pc -3 , one derives a rate 
of one capture every 10 7 yr. The time a MBH would 
have to wait before tidally capturing a star is, assum- 
ing a typical relative velocity at infinity v 2 cl = 6<r 2 , 



longer than the dynamical friction timescale by a factor 
In A (to^/to,) 1 / 3 (w*/<t) 5 / 3 , and shorter than the capture 

timescale of a stellar hole by a factor to^ 3 > 800. Be- 
cause of the high velocity dispersion, the binaries formed in 
galaxy centers will be much closer than those formed in the 
cores of globular clusters. Indeed, close encounters leading 
to the tidal disruption of the main-sequence star (initial 
pericenter separation r p < r t or v ic \ > u* f=a 600 kms -1 ) 
will be very common. Even when r p > rt, tidal effects dur- 
ing each subsequent periastron passage will tend to circu- 
larize the orbit, and the high tidal luminosity will cause the 
star to expand and overflow its Roche lobe. Thus we can 
reasonably expect that tidally-captured stars will undergo 
some form of rapid mass loss (disruption) or mass trans- 
fer. If the debris are accreted with high radiative efficiency 
by the MBH, then this may naturally lead to an off-center 
('non-AGN') source with a peak X-ray luminosity close to 
the Eddington limit, 2 x 10 40 (m./150 M©) ergs -1 , shin- 
ing over a timescale of order 1.4 x 10 5 (to./150 M©) -1 yr. 
Non-nuclear compact X-ray sources with such high lumi- 
nosities have been observed, e.g., by ROSAT (Roberts & 
Warwick 2000). 

Finally, we draw attention to the fact that MBHs sinking 
towards the center may be captured by the supermassive 
black hole (SMBH) in the nucleus and removed from the 
cluster as they diffuse into low-angular momentum 'loss- 
cone' orbits. They would then be swallowed either by di- 
rect scattering into the SMBH or by gradual shrinking of 
their orbits due to gravitational wave emission. The case 
for stellar holes has been discussed by Sigurdsson & Rees 
(1997) and Miralda-Escude & Gould (2000). The rate of 
orbit decay through gravitational radiation is proportional 
to to,, and so is the characteristic strain amplitude of a 
wave emitted by a MBH in a given orbit around a SMBH of 
given mass (for to. <C A/ sm bh). The capture of a 150 M© 
hole by a central black hole with M sm bh ~ 10 6 M© would 
be detectable by the planned Laser Interferometer Space 
Antenna (LISA) out to high redshifts. 

We have benefited from informative discussions with A. 
Ferrara, X. Hernandez, M. Salvati, and S. Woosley. Sup- 
port for this work was provided by NASA through ATP 
grant NAG5-4236 and a B. Rossi Visiting Fellowship at 
the Observatory of Arcetri (P.M.), by the Royal Society 
(M.J.R.), and by the EC RTN network "The Physics of 
the Intergalactic Medium" . 
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